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We consider the simple inflationary model with peculiarity in the form of "plateau" in the inflaton 
potential. We use the formalism of coarse-grained field in order to describe the production of metric 
perturbations h of an arbitrary amplitude, and obtain non-Gaussian probability function for such 
metric perturbations. We associate the spatial regions having large perturbations /i ~ 1 with the 
regions going to primordial black holes after inflation. We show that in our model the non-linear 
effects can lead to overproduction of the primordial black holes. 
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I. INTRODUCTION 



Starting from pioneering works by Zel'dovitch and Novikov [1], and also by Hawking [2], the primordial black 
holes (hereafter PBH's) were subject of extensive ivestigations. The presence of PBH's may significantly influence on 
physical processes and effects in the Universe (such as nucleosynthesis, CMBR spectral distortions, or distortions of 
7-ray background radiation) due to Hawking effect [3], PBH's may be a component of dark matter (see e.g. [4], [5]). 
The formation of PBH's is determined by small scale, but large amplitude inhomogeneities in the Early Universe, and 
the processes of PBH's formation, evolution and decay link the physical conditions of Early Universe with conditions 
in the radiation-dominated epoch and present-day cosmology. Even the very absence of PBH's may significantly 
constraint the models of the beginning of cosmological evolution. 

Usually the processes of PBH's formation are associated with production of the scalar mode of perturbations during 
inflation (see e.g. [5-9]) or phase transitions in the Early Universe [10]. In this paper we are going to discuss the first 
possibility, which allows to use the powerful and well-elaborated theory of instability of the expanding Universe for 
analysis of conditions, under which PBH's can form. 

The theory of generation of adiabatic perturbations during inflation started from pioneering papers [11-13]. It was 
established that the RMS-amplitude of metric perturbations S^ms is connected with the parameters of inflationary 
theory by means of relation 

™^ " 2^ 101 ' 

where H is the Hubble parameter, </> is the velocity of the field, evolving during inflation. To get PBH's abundance 
in an observable amount, one should have Srms ^ 10~^ — 10~^ (see, e. g. [14]). On the other hand COBE CMBR 
data, as well as analysis of Large-Scale Structure formation constraint the amplitude of perturbations Srms ~ 10~^ 
at super-large scales. Therefore to get PBH's one should increase the amplitude of the perturbations by a factor 
10'^ — 10'* at small scales. Unfortunately this cannot be reached in the simplest inflationary models, since in these 
models 5rms logarithmically grows with increase of scale, and one should use nonstandard models having additional 
power at small scales to obtain significant PBH's amount. 

Recently, several models of such type were proposed. For instance, Carr and Lidsey [6] proposed toy model having 
blue type spectrum (the spectrum 5rms{k) oc A;", where k is the wavenumber, and a is the spectral index), and 
investigated the constraint on the spectral index a associated with possible PBH's formation in such model. Linde 
[15] has shown that blue type spectra can be naturally obtained in the two-field model of so-called hybrid inflation. 



Another type of model having a spike in the power spectrum at some scale kth was proposed by Ivanov, Naselsky 
and Novikov ([5], hereafter INN) . They considered one-field inflationary model with inflaton and assumed that 
the potential V{(j)) has a "plateau" region at some scale kth, and has a standard form (say, power-law form) outside 
the "plateau" region. The field (j) slows down in the "plateau" region giving increase of the spectrum of perturbations 
at the scale k^fi according to cqn. (1). One can adjust the parameters of "plateau" region to obtain the desired 
increase of the spectrum, and consequently the desired PBH's amount. Garcia-Bellido et all [8] and also Randall et 
all [9] considered more realistic two field models having a saddle point in two-dimensional form of potential V{(j), tp). 
Like the one-field model, the evolution of the system of fields slows down near the saddle point giving an increase of 
the spectrum power. Randall et all pointed out that such models solve several fine-tuning problems of the standard 
inflation, and therefore look very natural from the point of view of high energy physics. Garcia-Bellido et all carefully 
investigated the process of PBH's formation in such models (see also recent work by Yokoyama, [18]). 

If the primordial black holes are not super-large, they probably collapse during the radiation dominated epoch of 
the evolution of the Universe. This means that the amplitude /i* of the metric inhomogeneities inside the regions going 
to PBH's should be of order of unity to overcome the strong pressure forces during collapse of the perturbed region 
[14]. These large amplitude metric inhomogeneities are assumed to be generated during inflation as rare events in the 
random field of the metric perturbations. Since the amplitude of the inhomogeneities ft.* is rather large, the natural 
question appears: to what extent we can rely on the linear theory of perturbations which usually gives Gaussian 
probability distribution of PBH's formation? 

To answer this question we can apply the formalism of coarse-grained fields (introduced by Starobinsky [19]) as an 
alternative approach to the linear theory that can describe large amplitude deviations of the field and the metric from 
the background quantities. According to this approach, the spatially inhomogeneous field (j^ix, t) is divided into two 
parts: the large-scale part 0;^, which consists of the modes with physical wavelengths A oc ak~^ greater than some 
characteristic scale Xc-g > H^^, and the small scale part which consists of modes with A < Ac-g. During inflation, 
the physical wavelengths are stretched, and new perturbations are added to This effect may be considered as a 
new random force f{t) in the equation of motion of the field (fUs, and usually the dynamics of is described in terms 
of diffusion equation for probability density ^'((/>is, t). This equation was subject of a number of works in connection 
with problems of Quantum Gravity and Large-Scale Structure formation. Recently, it was pointed out, that this 
equation can be employed for calculations of the probability to find large amplitude peaks in the random distribution 
of field (pis, and it was mentioned that such approach can be applied to the problem of PBH's formation [20]. 

Here we would like to note that when studying the effects originating after the end of inflation, such as PBH's 
formation, one should use the large scale part of metric instead of large scale part of field. Contrary to the field 
(pis, the large scale part of the metric, namely the "inhomogeneous scale factor ais{x)" (see eqns. (23 — 24) for exact 
definition) is the quantity conserving during the evolution outside the horizon, and this property allows to connect 
the physical conditions during inflation with the physical conditions during radiation-dominated epoch, when PBH's 
are formed. Moreover, the criterion for PBH's formation can be directly formulated in terms of ais{x) ( Refs. [21], 
[22]). Therefore, the calculation of ais{x) gives a tool to describe quantitatively the generation of non-linear metric 
perturbations, and the evolution of these perturbations into PBH's. 

In this paper we calculate the probability distribution function V{ais{x)) in the model with almost flat region in 
the inflaton potential. The main idea of our calculations has already been applied in the models of so-called stochastic 
inflation (see, e. g. [23] and references therein), and is very simple. When the field (pia evolves inside the plateau 
region it slows down, and the random kicks (described by the force f{t)) significantly influence on its evolution. So, 
the trajectory of the field inside the plateau region becomes stochastic, and the time Ai that the field spends on the 
plateau, depends on the realization of the stochastic process. The total increase of the scale factor ais during the 
field evolution on the plateau, is obviously determined by At: ais oc e^^* . Since different regions of the Universe 
separated by distances greater than evolve independently, the increase of ais corresponding to different regions 
is determined by different realizations of the random process. Thus the scale factor ais varies from one region to 
another after the field passes the plateau, that is the quantum effects generate the coordinate dependence of the scale 
factor. The shape of function ais{x) is conserved during the subsequent evolution of the Universe until the scale of 
inhomogeneity crosses horizon at the second time. At that time, in the regions with significant contrast of ais (x) the 
primordial black holes are formed. 

Using the approach described above we calculate the probability distribution function V{ais{x)). With help of a 
simple criterion of PBH's formation we relate 'P{ais{x)) to the probability of PBH's formation. We show that in our 



^See also the papers by Hodges and Blumenthal, Hodges et all[lQ\ and Kates et all [17], who employed similar models in 
contest of Large-Scale Structure formation theory 
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case the non-linear effects over-produce PBH's q Although this result is very important qualitatively, it does not 
significantly change the estimate based on the linear theory. 

We use the simple one-field model, proposed by INN (see also Refs. [24], [25]). Due to simplicity of this model the 
bulk of our results are obtained analytically. We hope that our approach provides a reasonable approximation to the 
case of more complicated two-field models. We are going to discuss these models in our future work. 

The paper is organized as following. We introduce our model and discuss the classical dynamics of the metric and 
field in Section 2. In Section 3 we obtain an expression for V{ais{x)). We consider the role of non- linear effects on the 
statistics of PBH's production in Section 4. We summarize our conclusions and discuss applicability of our approach 
in Section 5. 



II. THE DYNAMICS OF CLASSICAL MODEL 

In this Section we consider the classical dynamics of spatially homogeneous parts of metric and field in the simplest 
inflationary model with a single scalar field (inflaton) and with the peculiarity in the inflaton potential. In this case 
the system of dynamical equations contains only two dynamical variables -scale factor a(t) and spatially homogeneous 
part (f)o{t) of the field 0, and reduces to the Hamiltonian constraint equation 



and to the equation of motion for field 



3i70o + = 0, 



where H — and another symbols have their usual meaning. We hereafter use the natural system of units. 

We assume that the effective potential V{(j)) has a small almost flat region ('plateau') between some characteristic 
values of field (pi and 02 • The potential is also assumed to be proportional to 0^ outside the 'plateau' region 

Vi^) = ^ 4 

at < 01 

^(0) = l/(0i) + A(0-0i) 5 

at 01 < < 02, and 

F(0) = ^ 6 

at > 02. Here l^(0i) = ^j^, A — A(^) -f- ^AL^h^J^^ As we will see below the size of the flat region is very small 

A0 = 02 — 01 <C 0, ^ 1 so we can set A « A. At sufhciently large values of 0o > 1 the kinetic term in the 

equation (2) is negligible in comparison with the potential term 

I « ^(0o), 7 
and the equation (2) reduces to an algebraical relation between H and 0o (so-called slow-roll approximation) 



From the equation (8) it follows that the Universe expands quasi-exponentially {H « const and a oc e^*) at 0o > 1 



^Note that this result differs from that obtained in Ref. [20]. 
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It can also be easily shown that outside the plateau region the field moves with large friction at > 1, so 

|.^o| « |3if^|. 9 

The friction dominated condition (9) helps to simplify the integration of the system (2 — 3). Integrating the eqns. 
(2 — 3) with help of inequalities (7), (9) at > </'2, we have 

(t>o {t) = ^0 exp - i^l^t) , 10 

and 

a{4>o) = ao exp {N{4>a) - A^(0o)), H 
where (po and ao are some initial values of the field and scale factor. 

N{ct>o)= r° Hdt = n{(t>l-4>1) 12 

is the number of e-folds of the scale factor during the field rolling down starting from some initial value of and down 
to the field (^2- The similar formulae hold at (l)end < (po < 



a{(f)o) = ai exp(iVe„d((/)i) - Ne„d{(l>)), 14 

where 0o(ii) = 4'i s-nd ai ~ and A^erid(0o) is the number of c-folds up to the end of inflation: Nend{4>a) — 

7r((?!)o — 4'lnd)^ where we assume that inflation ends at standard (for X(f)^ theory) value of (pend = Note, that 

Nendi4>i) should bc rather large. For example, to get a feature in the spectrum at scales, corresponding to the solar 
mass, we should have Nend{(f>i) ^ 50 — 60. Therefore, the value of cpi should be greater than unity {(pi 4.5 for 

Nendi^l) ~ 60). 

Now let us consider the dynamics of inflaton in the "plateau" region (pi < cpo < (p2- In this region the equation (3) 
is simplified to 

^0 + SHo^o + ^ = 0, 15 



where Hq = -^J ^Vq. The solution of eqn. (15) can be written as 



and for the field velocity we have 



kne-^"°' - 17 



where <pin = (/)oUo=02 = ~35^^^('^2) = 'i^'f'^ ^^'^ field velocity at the moment i=0 of entrance of the field in 
the "plateau" region. The second term in the eqn. (16) and the first term in the eqn. (17) are due to inertial influence 
of initial velocity and the last terms in the both equations are due to nonzero slope of potential in the plateau 
region. The evolution of the field in the plateau region can be divided into two stages. At first stage the field evolves 
mainly due to inertial term, and velocity exponentially decreases with time. After some characteristic time the 
nonzero slope of potential A starts to determine the evolution, the velocity tends to the constant value (j>fd = 
and the field amplitude starts to decrease linearly with time. The time can be estimated by equating the inertial 
and potential terms in the eqn. (16), and is determined by the condition 3Hot»,e^^°*' = ^, where B = ■^V{(po = (pi)- 

As we discussed in Introduction, the spectrum amplitude is inversely proportional to the field velocity {Srms ~ ^"^)) 
therefore we need to slow down the velocity approximately by ^ 10^ — 10* times to get the increase of the spectrum 
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amplitude from the initial value 5rms{in) = ^ 10~^ up to the typical for PBH production Smis ^ 10~^ — 10~^. 

For that, we should fix the "amplification" parameter a = ^ ~ 10^ — 10^. 

Our model has two possible limiting variants depending on the relation between the time tc of the crossing of 
plateau region by the field (j)o (00 (ic) = 0i) and t». If tc ~ the field crosses the plateau mainly due to inertia. 
In this case the parameter a determines the number of e-folds during plateau crossing 5N Hotc ~ ^Ina « 2.3, 
and therefore the width of produced bump in the spectrum remains small and fixed. The model of similar type was 
discussed by INN. Here we consider another possible case tc > , where the field spends some time on the plateau, 
evolving in the friction-dominated approximation. In this case the width of the spectrum is determined by the value 
of tc, which is the free parameter of our model. Instead of tc we will parameterize our model by the quantity 7-the 
ratio of wave numbers, corresponding to the fields 02, respectively, tc = Hq^ In 7. The parameter 7 cannot be too 
small 7 > a^^^ and we take 7 w 10'^ in the estimations. If 7 is not extremely large ln7 <C N{(j)i), the size of plateau 
A0O = 02 ~ 01 is of order of typical size Acj)^ = The typical relative size of plateau is very small 



97P 

A0O 



00 6^02 67V(0i) 



0.003. 18 



Thus, the correction due to the presence of plateau practically does not influence on the dynamics of the field outside 
plateau region and we can set A = A. On the other hand, the size of plateau is much greater than Hq - the typical 
size of quantum fluctuations, A0* — ^" ~ 10^ Hq- 

Typically, the estimate ^ 1 holds for arbitrary power-low potentials V{<p) cx 0^ provided power p is not very 
large. However the opposite limiting case is also possible. For example, Bullock and Primack [20] proposed the 
potential of the form 

y(0) = Af,p(l + arctan (0)), > 

F(0) = Abp(l + 4* 10^^021), 0<O 19 

where the constant Abp = 6 * 10^ is chosen to normalize the large-scale part of spectrum to the RMS-amplitude 
sa 3* lO^'"^. The flat region in this potential starts from = and ends at = —1.23* 10"^, and inflation ends itself at 
= (f'end = —1.55 * 10~2. It was mentioned by Bullock & Primack that this potential leads to strongly non-Gaussian 
statistics of field perturbations. 



III. NON-LINEAR METRIC PERTURBATIONS FROM THE QUANTUM DYNAMICS OF 

COARSE-GRAINED FIELD 

It is well known that there are two equivalent ways to describe inhomogeneous Universe. The first way is to consider 
inhomogeneities as a small corrections to the homogeneous space-time and study them in the frameworks of linear 
theory of perturbations. Another approach splits the metric and the field into large-scale part (coarse-grained over 
some scale greater than horizon scale), and small-scale part. During inflation, the dynamical equations for coarse- 
grained field 0;s and coarse-grained scale-factor ais are equivalent to eqns. (3, 8) provided the quantum effects are 
switched off. The quantum effects continuously produce new inhomogeneities of random amplitude with scales greater 
than the scale of coarse-graining. These inhomogeneities should be added to 0;s and ais and effectively this leads to 
the presence of stochastic force term in the equations of motion. Therefore, the dynamics of coarse-grained variables 
can be described in terms of the distribution functions of 0;^ and ais-, and in principal these distribution functions can 
provide the same information as the power spectrum of perturbations, and furthermore the coarse-grained formalism 
gives a tool for description of the metric perturbations with amplitude, greater than 1. 

The effective dynamical equation for the field 0;^ has the form [19] ^ 

0i, + 3Hu^u + ■i-V{<j)is) - D'^^fit), 20 



^See also recent papers [26] and references therein. 
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where D = . \% , and f(t) is delta-correlated random force, < f(ti)f(t2) >= S(ti — 12)- The equation for coarse- 
grained scale factor ais remains unchanged 



Hu = ^^V{4>u). 21 

The solution of the set of eqns. (20, 21) is extremely difficult problem, and can be done under some additional 
simplifying assumptions. For example if we choose the featureless potential, and consider the friction-dominated 
solutions of the cqn. (20), wc can obtain the solutions describing self- reproduced inflationary Universe (provided the 
stochastic term in (20) dominates over potential term, see for example Linde [23]). In our case we cannot use the 
friction-dominated condition in the beginning of the field evolution inside the plateau region. However we can adopt 
another simplifying assumptions: first we can set His = Hq = const inside and near the plateau region, and second, 
we can omit the stochastic term in the eqn. (20) outside the plateau region, assuming the field moves along the 
classical trajectory there. Under these assumptions the statistics of the scale factor ais is totally determined by the 
time that field (j)is spends in the plateau region 

A7V = ln{aoutlain) = HoM, 22 

where Um is the value of scale factor at the time t = oi entrance of the field in the plateau region, and Gout 
corresponds to the moment At, when the field leaves the plateau region. To sec; that let us consider the evolution of 
the scale factor ais in the comoving coordinate system. Outside the horizon the hypersurfaces of constant comoving 
time tcom practically coincide with hypersurfaces of constant energy density e = const. On the other hand, the field 
4>is evolves slowly during inflation and hypersurfaces of constant energy density arc close to hypersurfaces (fus = const, 
and therefore we can put ais{tcom) = (iisi4'is)- After the field passes the plateau region, the evolution of ais{4>is) can 
be described by the standard expression (14), so we have 

aisiM - flin exp {7T{<j>l - (1)1) + AN), 23 

where AA^ is nearly constant inside of coarse-grained regions with comoving scale Xc-g ~ UoutHQ^, but changes from 
one region to another. Thus, the metric outside horizon has the quasi- isotropic form 

ds'^ = dt^ - afsi4>o)aisix)Sjdxidx^ , 24 

where we represent the scale factor ais{(pis) as a multipHcation of two factors: a{(j)o) and ais{x) = e^^ . Here ais{(j)o) 
and (j)o{t) are determined by the classical equations (13), (14), and the spatial coordinates x are coarse-grained over 
the regions with scale \c-g- To estimate the change of metric from one region to another quantitatively, we introduce 
the definition of non-linear metric perturbation 

ais{(l)is) - a{(j)o) TTfA, 1 OK 

h = — = exp Ho{At - ic) - 1 25 

a{<Po) 

(remind, that tc = Hq^Ihj is the time which the field spends in the plateau region moving along the classical 
trajectory when the stochastic term in (20) is switched off). Note, that in the limit of small /i <C 1, the metric 
assumes the form 



= dt^ _ o^iMil + 2h{x))5]dxidx^ , 26 

and the definition (25) is reduced to the standard expression for growing mode of adiabatic perturbation outside the 
horizon. Namely, in this case h reduces to gauge independent quantities, introduced by a number of authors [11-13], 
[27] up to a constant factor. The variables (25, 26) do not depend on time outside the horizon. Therefore, using 
of these variables is very convenient to match the perturbations, generated during inflation with the perturbations, 
crossing horizon at the normal stage of the Universe evolution. As one can see from (25) the metric perturbations arc 
determined by stochastic variable At and the distribution of At must follow from the solution of eqn. (20). Note, that 
the definition of non-linear metric perturbations should be taken with a caution. In principal, one can use another 
definition relating to (25) by some non-linear transformation, and having the same limit (26) in the case of small 
h. For example. Bond and Salopck [28] used the quantity h = In ( ) to define non-linear metric perturbations. 

However, the criterion for PBH's formation can be directly expressed in terms of the quantity (25) (see next Section), 
and therefore this quantity is the most natural variable for our purposes. 



vi 



Although the assumption of constant Hq greatly simplifies the problem it still remains rather complicated for a 
simple analytical treatment ^ 

For further progress we have to make some additional assumptions. We will consider below the plateau region of 
sufficiently large size. For this case the field approaches to the end of plateau in the friction-dominated approxima- 
tion, which greatly simplifies the treatment of diffusion processes. To estimate the relevance of friction-dominated 
approximation we should compare the time tc and the time ^ In (a) of the decay of the inertial term (/) in the eqns. 
(15 — 17, 20). If tc > i* and therefore 7 ^ a^/'^, the inertial term in these equations can be neglected at < t < tc- 
In this regime the solution of the classical equation of motion (15) has the form 

0o(T)~02-aT, 27 

and the equation (20) becomes 

+ a = d^/^f{T), 28 



dr 



where /? = 3iJo, and we introduce the dimensionless time r — f3t, a = A/ (3'^ and d — — The stochastic 

equation (28) is associated with simple diffusion type equation, describing the evolution of positions probability 
distribution ^'(t, c^) 

a* _ ^ 92 ^ a ^ 29 
dr dcj) ' 

Now we assume that the distribution "if is not spread out sufhciently before t* — /3i* and take (5-distributed "if function 
at the moment r = as the initial condition for our problem 

*(rO = <5(0,, - (/.*), 30 

where (f>^ = Acj) — ar* is the value of field corresponding to the beginning of " friction-dominated" part of plateau 
region 0. 

Together with initial condition (30) we should specify the boundary condition at 4>is = (t>i ■ This condition depends 
on the form of the transition layer between the plateau region and the part of potential with steep slope -§^y{<t') = B. 
We assume this transition to be sharp, and therefore set the condition of absorbing wall at the downstream point 

4>is = 4>i 

^'(0i,r)=O, 31 

Note, that this boundary condition was used by Aryal and Vilenkin [32] for analysis of stochastic inflation in the 
theory with top-hat potentials. In that paper it was shown that the more reasonable smooth transitions between the 
flat and steep regions of the potential are unlikely to modify signiflcantly the resulting distribution. 
In our case the probabihty density P{t) of time r relates to the solution of eqn. (27) as 

V{t) = S\^,^^^,^d^^, 32 



where we deflne by 5* the probability current S — d-^'^ + . The conservation of the probability current allows 
to estimate the correction term to eqn. (32) due to nonzero 5'((/)i). Assuming that fleld moves along the classical 



trajectory after (pis = we have S{—4'i) ~ B / 0^'^ « 5(+0i) ~ d-g^'^ . Therefore the correction to the expression 
(32) is — a^^ ^ — 10^"* times smaller than the leading term. 



*In this case our problem is reduced to the first-passage problem for one-dimensional Fokker-Plank (Kramers) equation, 
associated with eqn. (20) [29]. The general solution of this problem demands too much formalism [30], and is not considered 
here. Note, however that the simple asymptotic estimates are still possible in this case [30], [31]. 

^The estimates show that the characteristic width of 'I'(t,) is of order of H and much less than the size of "friction-dominated" 
region 0* — 0i. 
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The conditions (30,31) determine the solution of eqn. (29). This solution can be found by standard methods of 
the theory of diffusion equations (see, e. g. Ref. [33]), and in our case has the form 

11 2 

V^47rd(T-T,) 4d(r-T*) 

(l-exp{-— (^-0i)(^,-^i)}), 33 
a(r — r*) 

Substituting (33) to the equation (32) we find the exphcit expression for V{t) 

V{t) = ] (^^) exp {- ^ (01 - 0, + air - r.)f}. 34 

V (47r(i(T - T*) r- r* 4o!(r - r*) 

The expression for probabihty distribution of metric can be readily obtained from (34). Using the eqn. (22 — 25) to 
express the time r in terms of h, taking into account the eqn (27) and the definitions of a, d, and assuming a > 0, we 
obtain 

where 6pi = |^ = adrmsiin) is the standard metric amplitude calculated for the plateau parameters, and 

Nci=hij-T^/3, N,t=ln{l + h)+Nci 36 

are the numbers of e- folds for the classical path (f)o{t) and for a random path (pis{t), which start at 0* = <p{t*) and 
end at ^i. 

When the perturbations are small Ngt — Nd ^ h<^l, the distribution (35) has the standard Gaussian form 

r{h) = PG{h)= \ exp{--j^}, 37 

and in the opposite case of very large metric perturbations /i ^ 1 and Nst ~ In/i > Nd the distribution V(h) deviates 
sharply from the Gaussian law and has the power-law form 

V{h) oc 38 

As seen from eqns. (35 — 38), the non-Gaussian effects over-produce the metric perturbations of high amplitude 
in our model. To understand this fact, let us discuss the origin of non-Gaussian effects in our model. There are two 
sources for such effects. First, note that the "effective dispersion" alj^ = ^w-ZV^t in eqn. (35) depends itself on the 
value of the stochastic variable Nst- Qualitatively, it can be explained as follows. In the linear theory the dispersion 
= Sp[Nci is proportional to the time spent by the classical background field (po on the plateau. In non- linear 
theory the coarse-grained field 4)is{t) plays the role of background field, and therefore the distribution of the family 
of neighboring to (j) = (pis {t) paths should be described in terms of the probability distribution with dispersion j , 
which is proportional to the time spent by field (pis on the plateau. Second, the amplitude of large metric perturbations 
h depends on Ngi exponentially {h ~ e^'=), so order of magnitude increase of Nst leads to exponential increase of h. 
Obviously, these two effects increase the probability of large amplitude metric perturbations. 



IV. PROBABILITY OF BLACK HOLES FORMATION 



Although the distribution (35) provides very important information about the geometry of spatial part of metric 
outside horizon, it cannot be directly applied to the estimates of PBH's formation. Indeed, the distribution (35) is 
formed by the field inhomogeneities with wave- numbers k in the range (Afc = [kmin ~ CLinHo < k < kmax ~ o,outHo])- 
The process of PBH formation is determined mainly by the field modes with wave- numbers (Sk w fcb/i <C Afc), where 
kbh is the typical PBH wavenumber. The modes with k < kbh compose the large-scale background part of metric 
at the moment of PBH formation, and do not influence on the formation of PBH's significantly. The modes with 
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k > kbh lead to high-frequency modulation of the perturbation with k kbh, which is also unimportant, provided the 
mode with k kbu crosses the horizon second time at the radiation-dominated epoch. Therefore, in order to obtain 
the probability of PBH's formation, we should subtract the contribution of the large-scale and small-scale metric 
perturbations. 

In general it is very difficult to separate the perturbations of a given scale in the frameworks of non-linear approach. 
However, we can estimate the probability density of the perturbations, corresponding to the smallest scale k^h ~ 
ttoutHo 0. For that we simply put Nd = 1 in eqns. (35, 36), assuming that the random process starts when the mode 
with wavenumbcr fci = e^^aoutHo crosses horizon. This procedure automatically subtracts the large-scale contribution 
of modes with k < ki. The small-scale contribution is also absent due to our absorbing boundary condition. We have 

11 x'^ 
'P{h) = —j= r7:rexp{-— 2- — — ^} 39 

from the eqn. (35), where x — ln(l -f ft,), and in the limit of small h we obtain again the Gaussian distribution 
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The distribution (39) has nonzero first momentum Mi = dhhV{h) — ^Spi (the lower limit of integration should be 
-1, since the metric perturbations with h < —1 are cut off). The contribution of Mi should be added to the background 
part of metric, and further we will use the renormalized metric perturbation hr — h^^Spi instead of h. The probability 
to find the metric perturbations with amplitude greater than some threshold value /i* P{h^) — dhV{h) can be 
estimated as 

J_. 2Spi{x. + l)'/' ^ 

^*' V2^^ x4x.+2) ^'^P^ 2Sl,{x^ + iy' 

where = ln(l + |(5p/ + /i*), and we assume ft* ^ 6pi. The same quantity, but calculated for the Gaussian distribution 
takes the well-known form 



PG(/i)«^^exp{--^}. 42 



V2^ft/''P^ 2<5p, 



The observed quantities (such as, e.g. the matter density of PBH's in different cosmological epochs) can be easily 
expressed in terms of the probability P{h^), provided the mass of PBH's and some criterion for PBH's formation are 
fixed. In our case the criterion for PBH's formation should give the information about the threshold value ft*. Since 
this criterion plays very important role, let us discuss it in some details. First let us note that PBH's are formed from 
high amplitude peaks in the density distribution which are approximately spherically-symmetric (see e.g. Ref. [34]). 
It can also be easily shown that the maxima in the matter density correspond to the maxima in the function ais (x) . 
The form of ais{x) totally specifies the number of regions going to PBH's as well as dynamics of the collapsing regions. 
Therefore we formulate the criterion of PBH's formation in terms of conditions imposed on the function ais{x). 

The first criterion was formulated by Carr in his seminal paper [14]. It was shown that an over-dense region forms 
PBH if the density contrast at the horizon scale ^ lies approximately within the limits ^ < ^ < 1. The first part of 
this inequality tells that the over-dense region should stop expansion before the scale of the region crosses the sound 
horizon. The second part requires that the over-dense region does not collapse before crossing the causal horizon, 
and consequently the perturbation does not produce a closed world separated from the rest of the Universe. Then 
the criterion for PBH's formation was improved by Nadegin, Novikov and Polnarev [21] (hereafter NNP), and also by 
Biknell and Henriksen [22] with help of numerical computations. The initial condition used by NNP was chosen as a 
non-linear metric perturbation having the form of a part of the closed Friedman Universe matched with the spatially 
flat Universe through an intermediate layer of negative density perturbation. The conditions for PBH's formation 
depend on the size of this part (i.e. the amplitude of the perturbation), as well as on the size of the matching layer. 



In this connection, let us note that the black holes of smallest mass should give the major contribution to the present fraction 
of black holes, provided PBH's spectrum is flat (Carr, 1975 [14]). 
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The smaller matching layer is, the larger the pressure gradients needed to prevent collapse will be. Therefore, the 
amplitude of the perturbation forming PBH must be greater in the case of narrow intermediate layer. In terms of our 
function a{x) the NNP criterion reads 



K = —-l> 0.75-0.9 
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where is the value of a{x) at the maximum of the perturbation and a_ is the same quantity outside the perturbed 
region Q The first number on the right hand side of (43) corresponds to the matching layer of size comparable 
with size of the over-dense region, and the second number corresponds to the narrow matching layer. Assuming the 
matching layer to be sufficiently large we take /i* — 0.75 as a criterion of PBH's formation. 

Once the criterion is specified, we can link the desired PBH's abundance P{Mpi,ii) w P{h^,pi,ii) with the parameters 
of our model. For instance, consider the model having the matter density of PBH's equal to the critical one (the 
density parameter flpbh = !)• In this model we have [3], [6] 



/3(M) = 10-«( — ) . 
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Equating the expression (44) to the probability function (39), we have the equation determining the amplitude Sipi 
required for PBH's abundance (44) as a function of Mpbh 



PiKpbhjSipi) = PiMpbh), 



45 



and equating the expressions (42) and (44) we obtain the analogous equation for determining the reference amplitude 
d2pi when the non-Gaussian effects are switched off. The solution of these equations is given in Fig. 1. 
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FIG. 1. We plot the dependence of plateau parameter Spi on PBH's mass Mpbh assuming that the PBH's abundance is given 
by the eqn. (44). The solid line represents the solution of eqn. (45) (i.e we calculate 5pi taking into account the non-Gaussian 
effects in this case). The dashed line represents Spi calculated in the standard Gaussian theory. The PBH's masses lie in the 
range: lO'^^M© < Mpth < lO^Mpth. The PBH's of the mass W'^^Mq ~ lO^^g should be evaporated at the present time. 
Actually, the abundance of these PBH's is limited much stronger than is assumed in our calculations. 



One can see from this Fig. that the quantities 6ipi and S2pi increase with increasing of Mpbh and Sipi is always 
smaller than 62pi- It means that non-Gaussian effects over-produce PBH's in our model (at least when the simple 
criterion (43) is used), and the slope of potential can be steeper than that required in the Gaussian case. Typically, 



the ratio is about 1.5. Say, for the case of M, 



pbh 



Mq, we have Sipi{Mq) w 0.089 and S2pi{Mq) w 0.134. We 



plot the probability function 7'(/i) for 6ipi{Mq) — 0.089 in Fig. 2. 



''in the linear theory the density perturbation at horizon scale relates to the metric perturbation by ^ 
Therefore the estimate (43) is in agreement with Carr's result. 



\h (see, e.g. [27]). 



X 



0.20 



0.15 - 




0.05 - 



-15 -10 -5 5 

In (Mp,,/Me) 

FIG. 2. The dependence of probability density ■P(/i) on the metric amplitude h. The non-Gaussian curve (solid line) is 
calculated with help of eqn. (39) assuming PBH's abundance /3(Mo) « 10~®. That gives 5ipi{Mq) w 0.089. The dashed line 
is the reference Gaussian probability density calculated for the same abundance. For that curve we have 52pi{Mq) ^ 0.134. 
The dotted curve represents the Gaussian distribution taken with Sipi{Mq) « 0.089. This distribution strongly under-produces 
PBH's, and in this case we have /3 ^ 10"^''. 

In this Fig., we also plot the Gaussian probability function Vaih) for S2pi{Mq) = 0.134 (dashed line) and the same 
quantity for Sipi{Mq) = 0.089 (dotted line). Comparing the curves that correspond to the same PBH's abundance, 
we see that the non-Gaussian curve is flatter having larger values of V{h) at large h. The values of the Gaussian 
curve with the same plateau parameter 5ipi{Mq) is smaller by many orders of magnitude than the values of the 
non-Gaussian curve in the case of large h. 

Finally, let us note, that the non-Gaussian effects does not modify significantly the estimates based on the Gaussian 
theory. As we have seen, the ambiguity in the choice of the plateau slope due to these effects is about 1.5. This 
ambiguity seems to be less than the ambiguity in other parameters and can be obviously absorbed by a small change 
of the potential slope. 



V. DISCUSSION 

We have demonstrated that the non-Gaussian effects related to the dynamics of the coarse-grained field (inflaton) 
and to the evolution of the large-scale part of metric over-produce large-amplitude inhomogeneities of the metric 
compared to the prediction of the Gaussian (linear) theory of perturbations. We have derived an analytical expression 
for non-Gaussian probability distribution for non-linear metric perturbations, and estimated the influence of non-linear 
effects on the probability of primordial black holes formation. Wc used the simple single field inflationary model with 
peculiarity in form of the flat region in inflaton potential V(^), and power-law slope of the potential outside the 
peculiarity region. The key point of our approach is in the using of inhomogeneous coarse-grained metric function 
a{x) instead of the coarse-grained field (pis as a basic quantity. This allows to match the physical condition of 
production of inhomogeneities during inflation with the "observable" quantities. 

Our results can be considered as semi-qualitative only. The uncertainties come fromthe phenomenological character 
of our inflationary model as well as from the oversimplified treatment of the process of PBH's formation. The 
uncertainties related to the choice of parameters of inflationary model are mainly due to unknown form of the potential 
between the steep and flat regions, and also due to our friction-dominated assumption in the consideration of the 
stochastic process. These uncertainties can be eliminatedwith help of numerical simulations of stochastic process in 
more realistic models of inflation. The ambiguities concerning the criterion of PBH's formation are mainly due to the 
one-point treatment of this process. Actually, PBH formation is nonlocal, and dynamics of collapsing region depends 
strongly on the form of the spatial profile of the density perturbation (see e.g. Refs. [22], [35] for discussion of this 
point). The form of the spatial profile can be studied by means of n-point correlation functions of the coarse-grained 
metric and field. Unfortunately, the formalism of n-point correlation functions is still not elaborated (sec, however 
the Ref. [35] for the first discussion). Note, that probably the infiuence of the spatial profile of the collapsing region 
may be taken into account by a redefinition of the threshold value /i*, and this value might be effectively less. In this 
case the role of the non-linear effects would be damped. 

Finally we would like to note that the form of the distribution (35) does not depend explicitly on the specific 
parameters of our model. This allows to suppose that similar distributions can be obtained in more complicated 
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models, say, in two-field models proposed in Refs. [8], [9]. We are going to check this very interesting assumption in 
our future work. 
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